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We review the subject of perfect state transfer; how one designs the (fixed) interactions of a 
chain of spins so that a quantum state, initially inserted on one end of the chain, is perfectly 
transferred to the opposite end in a fixed time. The perfect state transfer systems are then used 
as a constructive tool to design Hamiltonian implementations of other primitive protocols such as 
entanglement generation and signal amplification in measurements, before showing that, in fact, 
universal quantum computation can be implemented in this way. 



1. INTRODUCTION 

The study of quantum mechanical systems is no longer 
restricted to the characterization of the properties of nat- 
urally occurring systems. Rather, with the advent of 
Quantum Information, the manipulation and engineer- 
ing of these systems to suit our purposes has come to the 
fore. One of the simplest ways to understand these ma- 
nipulations is to consider not universal Quantum Com- 
putation, but simpler sub-protocols. The sub-protocol 
of state transfer was originally proposed by Bose pQ, the 
motivation being that in a quantum computer based on, 
for instance, a solid state architecture, interactions are 
typically local, but we want to apply entangling gates 
between distant qubits. While a sequence of SWAP gates 
suffices to bring those distant qubits together, this is po- 
tentially prone to massive control errors, and it would 
be desirable to remove the need for such stringent con- 
trol. One promising alternative is to integrate a 'fly- 
ing qubit' such as a photon within the same system in 
order to transfer qubits quickly and easily. This, how- 
ever, requires the successful implementation of two in- 
formation processing realisations rather than just one, 
which is already a daunting task. Instead, the idea is 
to achieve state transfer by providing a pre-fabricated 
unit, made from the same solid state technology but with 
severely limited capabilities (thereby easing the fabrica- 
tion) which takes as input a state in one location, and 
outputs it at another. In principle, since it is not neces- 
sary to interact with the device except at the input and 
output, it can, to some extent, be more isolated from 
the environment and less susceptible to decoherence. Al- 
though this was the original motivation, whether these 
schemes are ever practically realized is largely irrelevant; 
they have already provided a huge degree of insight into 
engineering more complex quantum protocols. 

In this paper, we provide a review of how quantum sys- 
tems can be engineered precisely for the task of quantum 
state transfer with, in some sense, minimalistic proper- 
ties; there should be no interaction with the system ex- 
cept at initialization and read-out, and the system Hamil- 



tonian should remain fixed in time. We start, in Sec. [2] 
by proving necessary and sufficient conditions for perfect 
state transfer in a one dimensional system, and showing 
that this is equivalent to the protocol of entanglement 
distribution. These conditions provide an infinite family 
of solutions, and we discuss how to design systems with 
specific properties, and demonstrate that one particular 
solution is optimal with respect to a variety of param- 
eters. The understanding of this basic mechanism can 
then be applied to a variety of different systems, such 
as harmonic oscillators, and allow for the possibility of 
long-range couplings. 

In Sec. [3j we consider what happens if more than a 
single excitation is present, enabling us to generate en- 
tanglement, or to encode such that initialisation of the 
state of the system is unnecessary. In Sec. |4j we give a 
brief treatment of errors, although, if used as a construc- 
tive technique, this section is largely irrelevant. Sec. [5] 
describes the only known class of solutions for perfect 
state transfer in geometries beyond a one-dimensional 
chain such that the distances of transfer drops off no 
faster that a polynomial of the number of system qubits. 
Finally, in Sec. |6j we show how perfect state transfer 
schemes can be be modified in order to generate differ- 
ent types of multipartite entanglement, such as GHZ or 
W-states, followed by the much stronger result of how to 
design a Hamiltonian to perfectly implement a quantum 
computation without any external control. 



2. PERFECT TRANSFER IN THE SINGLE EXCITATION 
SUBSPACE 

The study of state transfer was initiated by Bose [TJ, 
who considered a ID chain of N qubits (open bound- 
ary conditions) coupled by the time-independent Hamil- 
tonian (X, Y and Z are the standard Pauli matrices) 

2V-1 

-^Bose = J J]] (X n X n+ i + Y n Y n+ i + Z n Z n+ i). 
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His protocol started with a state |0) 8Ar_1 , and simply 
allowed the Hamiltonian to evolve the system for some 
transfer time to, in the hope of creating the output state 
|q^® - ^ e ^ ne UI1 ] tnown s tate \ip) has transferred 
from one end of the chain to the other. For N — 2,3, 
transfer can be achieved perfectly. For larger N, it can 
be shown that perfect transfer is impossible [2]. It was 
subsequently observed that by modulating the couplings 
along the length of the chain, perfect state transfer can 
be achieved for all N IH E] , and it is this technique 
designing fixed Hamiltonians requiring no further inter- 
action that we review here, i.e. although a vast array 
of techniques for arbitrarily accurate state transfer exist 
with a variety of control paradigms, we restrict purely to 
the case of perfect transfer with no active control. While 
we will primarily restrict to a chain of qubits, many of 
the results can be generalized to larger local Hilbert space 
dimensions [6], including harmonic oscillators [7]. 

The solutions provided here for perfect state trans- 
fer represent an extremal point in the space of possi- 
ble solutions since they require a minimum of interac- 
tion with the system. As some level of control is intro- 
duced, perhaps through an encoding procedure (a Gaus- 
sian wavepacket [8] or across multiple chains [9j), or by 
some temporal modulation of the system Hamiltonian 
[TUHH], this can trade against the engineering require- 
ments imposed here, and it should come as no surprise 
that a plethora of protocols for high quality (possibly 
non-deterministic) state transfer protocols have arisen. 
A review of some of these can be found in [3] . 

While there are also extensions of these results to 
topologies other than one dimensional chains, such as 
H3HI8] . these are all subject to the significant crit- 
icisms that the transfer distance is no better than 
0(log N), and they require very high coordination num- 
ber. These facts suggest that they are unlikely to be 
realistic in a physical device, which should perhaps be 
restricted to localized couplings on a lattice of between 1 
and 3 spatial dimensions, and very cumbersome for ap- 
plication in other settings. For this reason, we will hence- 
forth restrict mainly to chains, and only mention briefly a 
single example a of a more general network (Sec.[5|. Nev- 
ertheless, it is important to note that there are poten- 
tially good reasons for studying these networks. Firstly, 
they can typically be constructed with fixed (i.e. not spa- 
tially varying) coupling strengths, which could reduce en- 
gineering requirements, albeit at the cost of other engi- 
neering requirements. Secondly, these more general net- 
works admit the possibility of being able to choose which 
of several output nodes a state should be routed to, ei- 
ther when the network is first constructed [14] or, poten- 
tially, during the transfer 13 . Finally, there is the motiva- 




FIG. 1 An array of qubits coupled by a nearest-neighbor in- 
teraction. 

tion from classical network theory, that some networks, 
such as those based on circulant graphs, are extremely 
robust to imperfections. These results, however, have 
not yet been transferred to the quantum regime. 

Thus, concentrating on the most direct method of 
transfer, we consider a chain of qubits coupled with the 
nearest-neighbor Hamiltonian 



N-l 



N 



H — § J n (X n X n+ i + Y n Y n+ i) — B n Z ni 

n—l n—1 

which is known as the XX Hamiltonian. The coefficients 
J n and B n are real numbers. The analysis of state trans- 
fer is vastly simplified by the observation that 

N 

71=1 

which means that H divides into a series of subspaces, 
characterized by the number of qubits that are in the 
|1) state. This is referred to as the number of excita- 
tions. It is not necessary to restrict to a Hamiltonian 
with this subspace form [Ml EI]; some work has been 
done on Hamiltonians such as 



N-l 



N 

H' = J n Z n Z n+ \ + B n X n , 



n=l 



where perfect state transfer can be achieved [2D], and 
many similar properties to those which we will demon- 
strate can be recovered [22] . We will see how these cases 
are also handled by the same formalism in Sec. |3.6[ 

Most of the state transfer properties can be derived 
by considering the subspaces containing no more than 
1 excitation and, in this case there is a large variety of 
other models such as Heisenberg Hamiltonians and har- 
monic oscillators to which these results also apply. We 



will explore some of these extensions in Sec. 2.7 



We shall initially assume that the system of spins is 
prepared in the state |0) = |0)® , and that an unknown 
state \ip) can be introduced onto spin 1 perfectly. The 
first of these assumptions will be relaxed in Sec. |3.3[ al- 
though the second one is necessary in a fixed Hamilto- 
nian setting. In order to relax that constraint, one could 
consider time- varying coupling strengths J\ and Jn-i 



a the only construction that we are aware of which gives perfect in any number of spatial dimensions with only a very mild re- 

transfer in a non-ID network with reasonable transfer distance laxation of assumptions, and hence has a large transfer distance 

b The routing problem can be solved on a regular lattice topology |19| . 
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PHH31 [M|- % expanding \ip) = a |0) + (3 |1), we observe 
that the component of the initial state alO)®^ is just 
required to evolve to the final state a\0)® N , but since 
|0)' 8Ar is an eigenstate of H, this happens automatically 
(up to a phase factor). Consequently, it suffices to con- 
sider the task |1) lO)®^ 1 -> {O)^' 1 11), which allows 
the restriction to the single excitation subspace of H. 
This is described by an N x N matrix, with the basis 
elements of this subspace being given by 



|n) = |0)^- 1 |l> |0) 



Using this basis, the Hamiltonian is represented by the 
matrix 



Hi 



fBx Ji 

Ji B 2 J-2 

J 2 B 3 



V 



Bn-i Jn-x 
Jn-i Bn ) 



where we have neglected a term — J2n=i B n \, as it evi- 
dently only affects the relative phase between the |0) and 
|1) components of after the transfer, and this can 
always be corrected by using a unitary rotation. 



1. Entanglement Distribution 

While we will be discussing the ability to transfer an 
unknown state from one end of a chain to the other, it 
is worth being aware that this is equivalent to the task 
of entanglement distribution Here, we imagine that 
two parties, Alice and Bob, want to share a maximally 
entangled state and they each control one spin of a cou- 
pled spin network. 

Lemma 1. The protocol of entanglement distribution, 
wherein Alice prepares a maximally entangled state, and 
puts half of it on the spin network, so that Bob can re- 
move it at some later time to, such that they perfectly 
share a maximally entangled state, succeeds if and only if 
the spin network is a perfect state transfer network from 
Alice to Bob. 

Proof. Let us assume that Alice controls two spins, and 
1, where 1 is the first spin of the network. Bob controls 
spin B on the network. Evidently, if Alice prepares a 
state 



(|oo> 0)1 + |ii} 0jl )|o> 8iV - 1 /^ 



(1) 



and waits a time to, then if the network is a perfect state 
transfer network, it evolves to 

(|00) ^ + |11) ^)|0)^ 1 /V2, 

and the protocol is successful. On the other hand, con- 
sider an arbitrary network, acting on the state (|00) 1 + 



|11) x ) \$in) A/2- If the output state \^f) is to be maxi- 
mally entangled between spins and B, it must satisfy 



U ®Ub\V) = (|00) 



n: 



O.B, 



for some single-qubit unitaries Uq and Ub- However, 
since (U <g) Z7*)(|00) + |11)) = (|00) + 1 11) ) , we can take 
Uq = 1 without loss of generality. Therefore, |0) 1 must 
have transferred to U' B |0), and must have transferred 
to U B |1). Thus, by linearity, we conclude that the sys- 
tem would have to transfer any initial state (a \ 0) 1 + 

i9 |l) x ) |$ in > to the state U B (a \0) B +/3\1) B ) |$ out ) i.e. up 
to a correctable local unitary on Bob's spin, the network 
must perform as a perfect transfer network, at least for 
some initial configuration of the system. □ 



2. The Symmetry Matching Condition 

It is now our task to consider what properties the {Ji} 
must satisfy in order to allow perfect transfer from |1) — > 

\m- 

Lemma 2. In order to transfer an excitation from 
|1) ~* \N) i- n a IB N-qubit nearest-neighbor XX -coupled 
system with open boundary conditions, the Hamiltonian 



must be mirror symmetric, 
BN + i_ n for all n. 



V - 



72 



id B„ = 



Proof. Let the eigenvalues and eigenvectors of Hi be A„ 
and |A„) respectively. The initial and final target states 
are given in terms of the eigenvectors by 



N 



|1) = 5>„|A, 

n=l 
JV 

\N) = Y,Pn\K 



(2) 
(3) 



Since the target is 



,-iH lto |^ = e t4> |^ 



for some phase <fi and time to, it must be true that 

e" a "S„ = e*Pn Vn. 

In particular, this reveals that \a n \ 2 = |/3 n | 2 , and by 
raising H\ to an integer power, m, we can relate 



N 



(i| ht li) = J2 K 



{N\ H™ \N) 



n=l 



For m = I, this gives that B\ = Bn. For m = 2, we 
find that B\ + J\ = B% + J%_ x , and thus ,]\ = J%_ v 
Each time that m is increased by 1, one new variable is 
introduced on each side of the equation (either a B n or 
J 2 depending on the parity of m), and one finds that they 
must be equal, giving the required symmetry properties. 

□ 
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This proof can be extended to more general coupling 
terms of the form 

Re( J n )(X n X n+ i+Y n Y n+ i)+lm( J n ){X n Y n+ i — Y n X n+ i), 

and one finds the condition \J n \ 2 — |JAr- n | 2 . The effect 
of complex coupling coefficients is just to alter the ulti- 
mate arrival phase of the |1) component of \ip) relative to 
the |0) [25], so it suffices to restrict to considering real, 
positive, J n . 

Given that B n = Bpj + i_ n and J n = Jn-u, the Hamil- 
tonian Hi commutes with the symmetry operator 

N 

5 = 2 |n) ( N + l-n \ . 

n=l 

Thus, Hi further subdivides into symmetric and anti- 
symmetric subspaces, with eigenvectors |A* ) and |A°) re- 
spectively. 

Lemma 3. A necessary and sufficient condition for per- 
fect state transfer in a symmetric chain is that there 
should exist a time to and a phase </> such that e~ lt ° x ™ — 
e^ and e~ itoX " = — for all eigenvectors (A ra |l) ^ 

Proof. Consider the decomposition of |1) in terms of 
eigenvectors, and evolve it with the Hamiltonian. 




This immediately reveals the required conditions for any 
a n 7^ 0. It follows that all initial states \n) transfer to 
their mirror opposite S \n) = \ N + 1 — n ). Furthermore, 
if state transfer occurs in time to, then by symmetry there 
is a revival of the state on the input spin at times 2nto 
for integer n and perfect state transfer at all times (2n + 
l)to. □ 

This represents a remarkably simple condition on the 
eigenvalues of Hi ; one which is readily tested, and widely 
applicable (with slight modification, it need not be lim- 
ited to a nearest-neighbor Hamiltonian, nor to an exci- 
tation preserving Hamiltonian). Within the single exci- 
tation subspace of a nearest-neighbor model the a n ^ 
condition can be relaxed by observing that this never 
happens. For an eigenvector |A„) = ^2 m ^n,m\lR ), the 
A n , m obey the recursion relation 

^m^n,m+l (An ^n)^n,m J71— l^n,m— lj 

so if A n . m _i = \ n ,m = 0, then \ n . m +i — 0. The special 
case of m = 1 shows that if A n ,i = a„ = 0, then A n ,2 = 0, 
and so it follows by induction that the entire eigenvector 
is i.e. there is no eigenvector with a n = 0. Alterna- 
tively, we would find that somewhere in the chain, there 
is a J m — 0, but this automatically forbids state trans- 
fer. The only remaining problem is how to identify which 
eigenvalues belong to the set {A^}, and which to {A°}. 



Lemma 4. For tridiagonal matrices, with negative off- 
diagonal matrix elements and eigenvalues A„ (n = 
1 . . . N ), the number of sign changes in the eigenvectors 
| A n ) is n— 1. The eigenvectors are assumed to be ordered 
such that A„ < A„ + i. 

This is proven in [27] . As a consequence, for symmet- 
ric tridiagonal matrices, the symmetry of the eigenvec- 
tors alternates, as it is determined by the number of sign 
changes (an odd number means that the eigenvector is 
antisymmetric). Our matrix H\ is a tridiagonal matrix, 
but with positive coefficients J. As such, the ordering 
is reversed, with the maximum eigenvector being the one 
with no sign changes. Consider {A„} to be an ordered set 
of eigenvalues, A„ < A„ + i. The necessary and sufficient 
conditions for state transfer in a symmetric chain then 
become that 

A„ - A„_i = (2m„ + l)7r/t , (4) 

where to is the state transfer time, and m n is a positive 
integer (which can vary with n). This means that the ra- 
tio of pairs of differences of eigenvalues must be rational. 



3. Inverse Eigenvalue Problems 

In Sec. |2.2[ we saw that in order to test if perfect state 
transfer can be achieved in a ID nearest-neighbor coupled 
chain, one only needs to test a symmetry condition, and 
then consider the eigenvalues of an N x N matrix. Rather 
than simply test such a condition, we would like to know 
how to fix the coupling strengths { J„} and {B n }. Fortu- 
nately, this is an extremely well studied problem. We can 
select any set of (non-degenerate) eigenvalues {A„} com- 
patible with Eqn. Q, and solve an Inverse Eigenvalue 
Problem |27j to find, in time poly (TV), the corresponding 
coupling strengths. Every possible choice of eigenvalues 
has a unique solution with positive (non-zero), symmet- 
ric, J n . 

Specifying a chain by its eigenvalues is mathematically 
elegant, but why should one select a particular set of 
eigenvalues? In the following subsections, we shall con- 
sider a specific set of eigenvalues, and show that they 
yield a family of chains which are optimal for a range of 
properties, such as transfer speed. For the remainder of 
this subsection, however, we shall outline a technique of 
Karbach and Stolze [2S\ that motivates how to select the 
eigenvalues so that the final chain has specific properties. 

The specific target that we have in mind is to create 
state transfer chains with as little spatial variation of the 
couplings as possible. One might hope that this would 
reduce manufacturing requirements for such chains . As 



c In practice, we will end up trading the degree of spatial variation 
for the required accuracy of the eigenvalues, and thus coupling 
strengths. 
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such, we start from a situation where all the coupling 
strengths are equal, J„ = 1, and B n = 0, which would 
be our ideal situation, except that perfect transfer is im- 
possible. The eigenvalues are readily found: 



A, 



-2 cos 



N + 1 



Let us denote min„(A„ — A„_i), the minimum spacing in 
the spectrum, by 8. We now select a quantity 7r/< -C 8, 
which will determine the state transfer time. In the 
present case, this means that to ~ TV 2 . Each of the 
eigenvalues can then be changed by no more than ir/to 
such that the perfect state transfer condition of Eqn. Q 
is satisfied. Subsequently, an Inverse Eigenvalue Prob- 
lem can be solved to find the coupling strengths of a 
chain with that spectrum. Since the eigenvalues were 
only slightly perturbed, the coupling strengths are only 
slightly perturbed, and thus the resultant chain is close 
to being uniformly coupled, and close to having no mag- 
netic field. One can improve on this by imposing on the 
choice of spectrum that A n = — Aat + i^„, so the magnetic 
fields necessarily come out to zero. 

Lemma 5. By selecting the eigenvalues of a symmetric 
tridiagonal matrix to satisfy A„ = — \n+i-u for all n, 
the result of the inverse eigenvalue problem automatically 
satisfies B m = for all m. 

Proof. We will prove this by giving the essence of the 
inverse eigenvalue algorithm 29J d . Consider the function 



B'(X n ) = j~| (A« — A m ), 

m—l^n 



which is the derivative of the function B(X) — Yim=i(^~ 
A m ), the characteristic polynomial of the system, eval- 
uated at A„. This can be expanded in the case of 
A„ = — Ayv+i-n such that 

/ N/2 

B'(X n ) = — 2X n I Y\ (A m — A n )(— X m — A„) 

I N/2 

B'(X N+1 _ n ) = 2A„ I ]Q (A m + X n )(— X m + A„) 

\ m— l^n 

for even N. Odd N follows similarly, so that, for any N, 
B'(X n ) = (-l) JV - 1 S / (A A r+i-„). In it is proven that 
for a symmetric tridiagonal matrix, 



1 



(5) 



d While containing the essence of the solution to the inverse eigen- 
value, this particular method is not the most numerically stable. 



Taking the ratio of these terms, 

K| 2 _ (-l^B'iX 



N+l-n) 



|ajV+l-n| 2 

This imposes that 



B'(X n ) 



= 1. 



JV 



(1| Hl m+1 



71=1 



2m+l 



a,, 



for every integer m. We must now compare this to an 
expansion of (1| iJ 2m+1 i n terms of J n and B n . Every 
term in this sum must contain an odd number of B n with 
n < to, and only one term J 2 . . . J^ n _ 1 B m contains B m . 
Provided J n ^ 0, we can therefore solve iteratively to 
find that all B n = 0. □ 



4. An Analytic Solution 

Rather than relying on numerical solutions, it is often 
beneficial to have special cases enumerated. A number of 
coupling schemes with integer spectra have been found 
[H 1301 131] j although we will only concentrate on 
the simplest of these, which is based on the J x rotation 
matrix of a spin ^(N — 1) particle. This can be written 
as a tridiagonal matrix 

JV 

J X =IJ2 Vn(N-n)(\n) (n + l| + \ n + 1) (n|), (6) 



which is exactly the form of H\ (J„ = ^^n(N — n), 
B„ = 0). You may recall that the eigenvectors must be 
the \M X ) states, with eigenvalues M x = ±(N-1),^(N- 
1) — 1 ... — |(iV — 1). Thus, the eigenvalues all satisfy 
A n — A n _i = 1, and since [J X ,S] = 0, the conditions 
for perfect transfer are satisfied, with to = w. To see 
this in a more natural way, consider the evolution of a 
Hamiltonian 

AT-l 



acting on the initial state |0) 



0N-1 



it, the state evolves to |1 
states 



Clearly, in a time 
Now let us introduce the 



1 



E 



V(V) ^Xn- 1 



where w x is the Hamming weight of the bit string x. One 
simply has to verify that 

H \lp n ) = J n _x IV'n-l) + Jn IVVrl) 

to see that it takes on the claimed diagonal form of J x 
on the subspace of {|^n}}, and that the claimed perfect 
transfer is achieved. Alternative proofs can be found in 
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5. Optimality of Solutions 

Evidently, arbitrarily fast transfer can be achieved by 
just continuing to scale up the coupling strengths. A 
transformation J n — > nJ n for all n yields t — > t a /n. 
However, in practical settings, we are typically con- 
strained by properties such as the maximum coupling 
strength, or possibly the maximum eigenvalue. In such 
situations, which is the best set of eigenvalues to choose? 
We will now show that the optimal choice with respect 
to a number of conditions is that of Eqn. (jfjj, as proved 
in [331 El]. 

1. Transfer Time 

We would like to optimize the transfer time t subject 
to some constraint. The most easily applied constraint 
is that of the maximum eigenvalue. Alternatively, this 
can be viewed as, for a fixed to, how can we minimize 
the maximum eigenvalue? The minimum spacing be- 
tween eigenvalues is 7r/to, and this is achieved by the 
analytic solution previously given. Thus, the difference 
A max — A m j n is a minimum. While bounding this range is 
useful, it is not the most physical parameter to optimize 
with respect to. Rather, it would be preferable to find, 
for a given maximum coupling strength J maX ) what the 
minimum transfer time is. This problem was first con- 
sidered in [M]. We will concentrate on the case of even 
N, as the odd N case is somewhat more involved. 

We define Sk — Xk+i — Afc to be the difference in eigen- 
values, and 5 = min^f^. Note that 6 > it /to, and we 
wish to consider all sets of models with the same to. The 
Hamiltonian divides into symmetric and antisymmetric 
subspaces, with the eigenvalues satisfying 

Y,(K-K) = Tr(SH 1 ) = 2J N/2 . 
For any particular coupling scheme, one has 

Jmax > J N /2 >\NS> \N^/t Q . 

Moreover, for the analytic scheme, these inequalities 
completely collapse, J max = \N /to, and thus, for a fixed 
to, this represents the scheme with the smallest J max . 

2. Timing Errors 

One of the underlying assumptions to the process of 
state transfer has been that at t = 0, a quantum state 
can be placed on spin 1, and at t = to, it can be removed 
from spin N, instantaneously. In practice, this will not 
be the case, so it would be useful to know how to design 
the chain so that for a given tolerable error e, the possible 
time window for implementing the transfer of the state 
is as long as possible. Critical to this consideration are 
the quantities 

ln {t) = {n\e- iH ^\V), (7) 




20 40 Time 60 80 100 



FIG. 2 A comparison of the transfer fidelity (i/>in| Pout l^in) of 
the uniformly coupled chain of 31 qubits (blue), the almost 
uniformly coupled chain giving perfect state transfer (red), 
and the analytic solution (green). 



and, more specifically, 71 (t). Given the symmetry of 
the chain, optimizing this function will also optimize the 
equivalent function for the arrival of the state at the other 
end of the chain. Let us expand 71 (t) in terms of the 
eigenvectors {|A n » of H x , 7l (i) = £^ =1 {a^e^K Of 
course, 7(0) = 1, and this must be a maximum. Thus, 
how quickly it changes is governed by 

J2 n 

d -± =-£KlX = -0|fl?li>, 

ai *=o „=1 

and minimizing the rate of departure of the excitation 
equates to minimizing J\ + B\ for a given to- It thus 
seems likely that our analytic solution, where B\ = and 
J\ is only 1 / ViV of the maximum coupling strength, will 
be extremely robust. For example, Fig. [2] compares this 
chain with another following the technique of Sec. |2.3| 
Further confirmation arises from a closer examination of 
the expression J2n=i l a »| 2 ^n- The analytic solution has 
the smallest possible maximum eigenvalue. Furthermore, 
upon solving for the eigenvectors [30], one finds that the 
a n corresponding to the largest |A„| are exponentially 
suppressed in N, and only those with essentially no con- 
tribution to A 2 have a contribution in |a„| 2 , which thus 
shows that the analytic solution must be very close to 
optimal, given that none of the a n = 0. 

6. Beyond Nearest-Neighbour Couplings 

There is no special reason to restrict to nearest- 
neighbor couplings in the Hamiltonian H; in practice 
there are likely to be residual couplings with strengths 
that decrease with distance. The introduction of 
XX + YY terms between non- nearest-neighbors does not 
change the fact that the Hamiltonian decomposes into 
distinct excitation subspaces, and the symmetry match- 
ing condition of Lemma [3] still holds. The Inverse Eigen- 
value Problem for this wider class of cases is less well 
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studied, and the existence of solutions is no longer guar- 
anteed - the alternating symmetry property of eigenvec- 
tors is not guaranteed to hold, which leads to problems 
with selecting the correct eigenvalues. Nevertheless, by 



using the technique of Sec. 2.3 to start from a solution 
where eigenvalues can be calculated, a near-by solution 
can be expected to exist, and then an efficient algorithm 
has been formulated to find the coupling parameters |33) , 
as will now be described. 

We are given a Hamiltonian H(r) which satisfies 



JY 



= o, 



and is represented by a symmetric N x N matrix Hi (r) 
in the first excitation subspace. This Hamiltonian de- 
pends on N parameters {r,-}, which can be thought of, 
for example, as the positions of the spins when they are 
coupled by a distance-dependent coupling strength. The 
desired eigenvalues are contained in the N x N diagonal 
matrix A c . 

The algorithm starts with a first estimate rim for f. 
The matrix H±(r^) is diagonalized by Uq, 



■eA )U, 



where Ao is a diagonal matrix which encapsulates the 
errors in the energies, and e is a small parameter. The 
second step of the algorithm should use a vector rm = 

r^o) + s Sr, and it is our task to find the best choice for 

Sr. A unitary U\ diagonalizes the Hamiltonian, 



H 1 {F^) = U 1 A(l + eA 1 )Ul 



(8) 



Provided e is small, Uq and U\ should be close to each 
other. As such, U\ can be parametrized in terms of e, 



Ui = U Q (1 + ieQ)(t - ieQY 



(9) 



where Q is a Hermitian matrix containing information on 
the change in eigenvectors. This parametrization ensures 
that UqIIi — > 1 as e — > 0, and the unitarity of U\ can be 
seen by rewriting 

isQ = (1 + U\U )- 1 {1 - U\U Q ). 

Given that the eigenvalues of U\Uq must be of the form 
e~ 10 , the eigenvalues of Q are real (tan(0/2)) i.e. Q is 
Hermitian. Conversely, for any Hermitian operator Q, 
the eigenvalues can be identified with tan(0/2), and thus 
U\Uq must be unitary. The parametrization of Eqn. (JoJ) 
can now be substituted into Eqn. (|8|, and expanded in 



Note that one must select a consistent scheme for ordering the 
eigenvalues on the diagonal. 



terms of e. The terms for e° cancel, and the terms for e 1 
give 



^— ' or. 



U = AAi - AA + 2i(QA - AQ). 



The aim of the iteration should be to choose Sr such 
that Ai = 0. Since the diagonal elements of the final 
term, QA — AQ, are zero, all the eigenvalue information 
is encapsulated by the diagonal elements of the equations, 
while changing eigenvectors only affects the off-diagonal 
elements. Thus, the previous equation can be rewritten 
for just the diagonal elements, 



M.5r 



(10) 



where b is a vector of the diagonal elements of — AAo, and 
the n th column of the matrix M is given by the diagonal 
elements of 



t dHi 
dr n 



U . 



The solution to Eqn. (10) is the vector Sr, which gives 

the correct eigenvalues to 0(e 2 ). Provided Sr is small in 
comparison to tvq), this protocol iterates, squaring the 
error at each step. Hence, to achieve an accuracy of Eq, 
only O(log(eo)) iterations are needed. Since there are ef- 
ficient algorithms for solving Eqn. (10), and because the 



matrices to be diagonalized are symmetric (hence there 
are efficient diagonalization procedures, such as House- 
holder reductions the cost of each iteration scales 
polynomially, 0(N 3 ), with the number of qubits in the 
chain [35] , and the required parameters can be found with 
an efficient classical computation. 



7. Alternative State Transfer Systems 

In this section, we have derived the necessary and suf- 
ficient conditions for perfect state transfer on a chain of 
spins coupled by an XX Hamiltonian, and found an ana- 
lytic solution which is optimal with respect to properties 
such as state transfer time. These same conditions can 
be applied to a number of other systems which are exci- 
tation preserving, although it is important to note that 
in the next section, when we move to using many excita- 
tions, the equivalences no longer hold. 

The first case that we shall consider is a modulated 
anisotropic Heisenberg Hamiltonian 



JV-i 



H\ 1C i S = | J n (X n X n+ i +Y n Y n+ i + A n Z n Z n+ i) 

n=l 
N 
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Within the first excitation subspace, the original matrix 
Hi can be recovered by setting 

b n = B n + |(J„A„ + J„_iA„_i), 

up to an identity matrix. One consequence of this, how- 
ever, is that it is impossible to achieve perfect state trans- 
fer with all the b n set to 0, unlike lemma 5] [55] . 

The other case of interest is a chain of harmonic oscil- 
lators [7], although this can be further generalized to a 
chain of g-deformed oscillators [37]. As before, there is a 
set of N sites, and at each site n we introduce the creation 
and annihilation operators a n and a n , which obey the 



bosonic commutation properties [a Tl 



, a. 



ISr, 



and 



[a n ,a m ] = [ajjjdtj = 0. A nearest-neighbor Hamiltonian 
of the form 



N 



N-l 



Hho = 2 ^2 B n a) n a n + ^ <^( a n+i a « + a n a n+l) i 11 ) 



is also excitation preserving, 



N 

Hho, 



= 0. 



Starting from the vacuum state, |0), the single excitation 
subspace can be defined by |n) = a n |0), and the Hamilto- 
nian within this subspace is the same form as H\. Thus, 
all properties of perfect state transfer arc recovered. 



3. HIGHER EXCITATION SUBSPACES 
1. The Jordan-Wigner Transformation 

The primary benefit of choosing to examine the 
nearest-neighbor XX model rather than the Heisenberg 
model is that when we come to analyze the presence of 
many excitations on the chain, vast simplifications arise. 
This is because one can apply the Jordan-Wigner trans- 
formation 38J, which maps the excitations on the spins 
into non-interacting fermions. This mean that the indi- 
vidual fermions behave exactly as they did in the single 
excitation subspace, and we only have to take care of the 
exchange phase. 

The Jordan-Wigner transformation is as follows. We 
introduce creation operators a n for fermions at site n, 
and identify them with 



n 



Z m (X - iY), 



Thus, one can rewrite the Hamiltonian H as 

N-l N 

Hjw = ^2 4(4 tt «+i + a n+i a n) + 2 ^2 B n a n a n , 

n—l n—1 

up to an irrelevant 1 term. Let's now take the ground 
state of the system, |0), the state with no fermions. The 



eigenvectors of Hjyy in the one fermion sector can be 
written as 

N 

^2 A «,™ a L 12) 

m=l 

with eigenvalues A„. Our claim that the Hjw corre- 
sponds to non-interacting fermions needs to be shown 
by diagonalizing the full Hamiltonian, not just the one 
fermion subspace. We therefore introduce new operators 
b n , and propose that they satisfy 

N 



b n ^ ' 'Vi,m fll m> 



m— 1 

with the intent to show that 

n 

Given that the {X n , m } define the eigenvectors in the first 
excitation subspace, they must satisfy 

This helps when we expand the b n in J2 n ^nb n b n , 

N 

n,m,/c— 1 

Now we use an orthogonality relation for the eigenvec- 
tors, 

JV 

(0| a m a\ |0) = ^2 ^n,mKi,k = hn,k, 

which reveals that 



71=1 



jV 



E 



^•nb n b n ^ y B m Cl} rn (l m -f- Jra — 1^7n^7n— 1 ~t~ Jrn.Q \n ^777+1 



n m—1 

= Hj W . 

Therefore, the {b n } describe independent fermions. 
Heuristically, this means that if excitations on either 
spins 1 or 2 (say) transfer perfectly, then a pair of exci- 
tations on both spins 1 and 2 transfer perfectly to a pair 
of excitations on spins N — 1 and N because they are 
independent. However, because the two fermions have 
exchanged (the one on spin 1 has gone past the one on 
spin 2), a multiplicative factor of —1 is also introduced. 

The wedge product [33J ED is perhaps the most useful 
tool that we can bring to bear as a result of the Jordan- 
Wigner transformation. This is used to combine states of 
single excitations into valid states of multiple excitations. 
For example, |1) A |2) denotes excitations on spins 1 and 
2. It has two crucial properties that convey the fermionic 
nature of the excitations. The first is exchange, |1) A 
1 2) = — 1 2) A |1), and the second is the exclusion principle, 
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ID A |1) =0. As a result of the exchange property, we 
must pick a standard ordering for terms; the excitations 
are written in order from left to right, so a state of 5 
qubits 1 01001) along a chain is written as |2) A |5). 

From these basic properties of the wedge product, fur- 
ther properties are readily derived. Consider two vectors 

l a ) = E„ a n N an d \b) = J2 n & ™ 1^)' which are correctly 
normalized. 

\ a ) A \ b ) = ^2 {a m b n - a n b m ) \m) A \n) 

The normalization of the new state is 
\ }Ja m bn ~ a„b m \ 2 = ^(\a m \ 2 \b n \ 2 - a m b* m a* n b n ) 

= 1-| (b\a) | 2 . 

Again, we find that the two excitations cannot be in the 
same state, \a) A \a) = 0, although this is no longer re- 
stricted to the computational basis. The eigenvectors 
in, for example, the second excitation subspace, can be 
calculated from |A n ) A |A m ), which is equivalent to eval- 
uating the Slater determinant [30] . The corresponding 
eigenvalues are A n + A m . As such, the evolution due to 
the Hamiltonian H can be found from 

(e- lHlt \a)) A (e~ lHlt \b)) . 

To verify that this statement is consistent, we expand 

l a > = En a ™ l A ™> and \ b ) = En b n 

e- im {\a) A |6)) = e~ iHt ^ a n b m |A„) A |A m ) 

= ^a„6 m e-^ A " +A ™)|A„)A|A m ) 

n,m 

= ]T ( a » e_U "' l A »>) A (b m e- iXmt |A m » 

n,m 

= (e~* Hlt \a)) A (e- lHlt \b)) . 



2. Entanglement Generation 

One of the most trivial applications of multiple exci- 
tations on a spin chain is in the dynamic generation of 
entanglement. Consider starting a perfect state trans- 
fer chain in the state |+) |0)® 2 |+), where |+) = 
(|0) + |l))/v / 2, and allow it to evolve for a time to- In that 
time, the two |+) states exchange positions. However, 
only when two excitations are present do they generate 
an exchange phase. Thus, the system evolves from 

5 (IQ) + ID + \N) + ID A \N}) 

into the state 

|(|0) + |JV) + |1) + |JV)A|1». 



Restoring the normal ordering of the operators, and ne- 
glecting the central N — 2 spins, the final state is 

i(|00) + |01) + |10)-|ll», 

which is a maximally entangled state between the two 
most distant points. This simple observation led to the 
proposal that perfect state transfer chains be used to gen- 
erate ID cluster states in a highly parallel process [^0] . 

3. System Initialization 

The work on multiple excitations can also be used to 
remove the requirement of initializing the spin chain in 
the |0) state before the transfer. Let us encode the un- 
known state \tp) into the first two spins of the chain, as 
the state a |01) + £ jlO) [3T]- We shall assume that the 
rest of the chain is in some unknown state \x), where x 
is a string of N — 2 bits denoting the position of excita- 
tions on the chain. If the transfer scheme works for all 
possible states \x), then it also works for all possible su- 
perpositions and mixtures, p. We are only interested in 
transferring the state |^), and thus, there is no require- 
ment to preserve the initial state of the chain. 

At the start of the state transfer process, we have the 
state 

(a|l)+^|2»A|f>3... JV , 
and this transfers in time to to the state 

(a\N)+P\N^l))A\x T ) 1N _ 2 , 

where |a; T ) is the mirror of state \x), incorporating any 
phase of ( — 1) that may result. Reasserting the normal 
ordering of the wedge product, this gives 

(-I)"- \S T )^ N _ 2 A (a | A) + p \ N — 1) ), 

and therefore the encoded state \ip) is transferred per- 
fectly. For general initial states p on the rest of the 
chain, the output state is not the re-ordered form of p 
because of the different phases appearing on the vectors 
\x T ), but this is of no concern to us. Heuristically, the 
process that we have gone through is to note that the 
only problem when transferring states in higher excita- 
tion subspaces is that exchange phases are generated, and 
may cause entanglement. However, this entanglement is 
generated between different excitation subspaces (yield- 
ing different numbers of exchange phases), and thus, by 
encoding within a subspace of fixed parity of excitation 
number, the state to be transferred does not become en- 
tangled. 

Although encoding this state requires a (nearest- 
neighbor) entangling operation, the eventual read-out 
does not; one can measure spin N — 1 in the X basis, 
and remove the state from spin N (up to a Z rotation, 
depending on the measurement outcome). 
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In Sec. 6.3 we will describe an alternative transfer and 



scheme which does not require any encoding of the quan- 
tum state to be transmitted, and is nevertheless insensi- 
tive to the initial state of the chain. This is achieved by 
a unitary rotation applied to existing transfer schemes 
such that the unitary maps the encoded states described 
in this section to un-encoded states in the new scheme. 
The transverse Ising model 20 , 22J also achieves this feat 
of not requiring an encoding by virtue of a special feature 
of the way it can be mapped into the XX model that we 
considered here, as we will see in Sec. 3.6 (thereby ex- 
plaining the close link between the coupling strengths 
and magnetic fields found in (20] [55] and the analytic 
solution we gave in Sec. 2.4). On the other hand, the 
bosonic system of Eqn. |11| does not have an entangling 
operation between different particles and hence works di- 
rectly in higher excitation sub-spaces, without the need 
for an encoding. 



4. Transfer Rate 

The state transfer protocol involves placing an un- 
known state \ip) on spin 1, waiting a time to, and find- 
ing the state perfectly transferred to spin TV. With the 
added ability to consider higher excitation subspaces, we 
can now consider the state transfer chain as a channel 
for communicating quantum states [33J , and then the in- 
teresting question is what rate of transmission of single- 
qubit states can be achieved? 

Let us start with a state \ifji) = «i |0) + (3\ |1) on spin 
1, and let it evolve for a time t r . At this time, we will 
introduce a state \ijj 2 ) onto spin 1 (tracing out the state of 
that spin). After further evolution for to — t r , we might 
hope that a high quality version of would appear 
at spin N, essentially unaffected by the introduction of 
\ijj 2 ). If so, this would allow transmission at a rate l/t r , 
yielding a significant improvement over the original 1 /to 
rate. 

At time t r , before \if>2) is added, the system is in the 
state ai |0) + j3ie~ lHltr |1). The state can be split into 
parts which are |0) and |1) on spin 1, 



-iHxt r 



ID 



-iHxt r 



ID 



7l(*r)|l»+7l(*r)|l) 



which allows spin 1 to be traced out. The trace gives a 
mixture of two components, 



at |0)+/3i(e 



-iff it,. 



1) 



7l(*r)|I» 



and 



0171 (*r) ID- 



This immediately reveals that the protocol can only be 
perfect if 7i(£ r ) = 0. The state \ip 2 ) is subsequently 
added to spin 1, and the system is evolved for time to — t r . 
In that time, the two contributions become 



ai a 2 |0) + aajSi \N) + ( ai [3 2 - a 2 l3i-/i{t r ))e~ lH ^ to - t '- 
+/3i/32(e-^ l(t °- t - ) |D) A \N) 



ID 



(a 2 |0) + [3 2 e 



-iHx{t -t r ) 



I»7l(*r)/3l 



From this, one can calculate the reduced density matrix 
p of qubit N, and then evaluate the fidelity 

F = (i> 1 \p\i> 1 ), 

The full expression is too complex to reproduce here, 
but one can more readily understand an example case, of 
Pi = 02 = 1, where it turns out that 

F=l-| 7l (^)| 2 -|7i(M| 4 - 

For the analytic solution of the state transfer chain, "fi(t) 
is a decreasing function of time (0 < t < to), with 



7i W 



cos 



JV-1 



irt 
2U> 



As a result, it is possible to achieve a transfer rate of 
\/N /to for some constant error probability e on each \ip n ) 
that is transferred. 

Since our focus is on the perfect realization of these 
protocols, we will focus on the case 71 (t r ) — 0. More- 
over, to achieve the target rate of l/t r , we will require 
that 7i(ni r ) = for all integers n = 1 . . . 2t Q /t r — 1. It 
is not always possible to enhance this rate; the smallest 
non-trivial case is N = 3 but, as a perfect state transfer- 
system, this is uniquely specified (up to an overall nu- 
merical factor of 1/to ) if we assume no magnetic fields 
can be used. One can readily show that 71 (t < to) > 0. 

Given that 71 (t r ) — implies that an initial state |1) is 
orthogonal to the state at time t r , the Margolus-Levitin 
Theorem [5T] can be used to bound the achievable rate. 
This theorem states that the minimum time for one state 
to evolve into an orthogonal one must be at least 



t r > 



Hi? 2 ID 



(l\H\l) 



7T 

2JD 



(12) 



thereby giving an upper bound to the rate of perfect state 
transfer for any given scheme. The rate is also trivially 
bounded by l/t r < N/2to, since it is impossible to store 
the states of more than N qubits on an A-qubit chain. 
The following lemma allows us to prove a tighter bound 
of l/t r < N/Ato, although we conjecture a stronger con- 
dition; that there are no chains with t r < to- 

Lemma 6. // a set of eigenvalues is chosen to fulfill the 
perfect state transfer condition of Eqn. Q), then a neces- 
sary and sufficient condition to perfectly achieve the rate 
M/2t for integer M is that all the R k for k = . . . M -1 
should be equal, where 



R k 



N 

E 



B'(X n ) 



(13) 
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and the sum is restricted to those terms satisfying the 
condition 



tc 



(A n — Ai) mod M = k. 



Proof. Using Eqn. (|5]), we can write that 



TV 



7iW = E 



-iX„t 



{ (-l)"B'(A n )£ 



/v 



= 1 (-l)™B>(\ m ) 



We will now demand that 7i(2mto/M) = for integer 
m = 1 ... M — 1 and 7i(2< ) = 1- Using Eqn. ([13]), the 
terms in the sum can be split up to give 



7i(2mt /M) 



EM— 1 p -il-KkmlM 
k=Q n k e 

l^k=0 H k 



Using the Fourier Transform, these equations can be in- 
verted to give that 



M-l 



R. , 



M E R ^ 



k=0 



as required. 



□ 



Note that if M is a composite number M1M2, then 
if the condition fails for either Mi or M2, it necessarily 
fails for M, since each of the R k for Mi is given by a 
sum of M 2 of the R k for M. However, the converse is 
not true. We conjecture that it is impossible to fulfill the 
condition of Lemma [6] for any M > 2, although we only 
have a proof for M > N/2. In the case of M = N, since 
none of the terms Rk can be if they are to be equal, they 
must all be R k = {-l) k /B'(\ k ). Given that B'(Ai) ^ 
—B'(\2) for iV > 2, it is immediate that it is impossible 
to perfectly achieve the rate N/(2to). Similarly, for any 
M > [N/2\ + 1, the different B'(X n ) must be partitioned 
between all R k , and it follows that at least two terms 
must be of the form (— l) n /B'(X n ), and hence different. 



5. Sequential Quantum Storage 

While we may have no way to construct a perfect 
state transfer chain with high rate (retaining the per- 
fect transfer property), it turns out that one can design 
a non-symmetric chain for which (1| e - iJI ^ n |1) = for 
n = 1 . . . 2to/t r — 1, which means that the states can be 
read in sequentially at spin 1, and read out sequentially 
from the same spin at a time 2to later. Recall that 



7i(«) 



N 

Ei 

71=1 



-iXnt 



In order to get periodic revivals at the input spin, then 
even without symmetry, the chain must have eigenval- 
ues such that the ratio of differences are rational (fol- 
lowing an identical argument to that of Lemma |3j). So, 



we can select these to be the simplest conceivable case, 
2t \ n /TT = -N + 1, -N + 3 . . . N - 3, N - 1. In the pre- 
vious subsection, we realized that if t r — 2to/N, then 
\ce,i\ 2 = l/N. To verify that this will work, we calculate 



71 (mt r ) 



1 

N 



N 

E< 



2N 



0. 



for any integer to which is not a multiple of N. With the 
{|a n | 2 } and {A„} specified, an iterative scheme for decid- 
ing the coupling strengths is evident; we simply calculate 

aitfrii} = EKi 2A ™ 



for all integers to, and solve for the coupling strengths. 
It can be verified that a matrix with B n = and 

2 _ n 2 (N -n)(N + n) 
n ~ (2n- l)(2n + l) 

has the desired properties, and N single-qubit states can 
be sequentially stored in the chain, which must be the 
optimal solution in terms of density of storage. 

Provided the bits are read-out in the reverse se- 
quence to which they were stored, there are no unde- 
sired side-effects from the exchange of fermions, because 
no fermions are exchanged. On the other hand, if one 
reads them out in the same sequence as they were read- 
in, there is a controlled-phase gate enacted between ev- 
ery single qubit. This could be a very useful feature for 
generating GHZ states - it would suffice to set all the 
states \ip n ) = |+), and the set of states that are out- 
put are equal, up to local unitaries, to an A^-qubit GHZ 
state. By being selective over the order of removal, one 
can control which spins get entangled with which others' . 
For storage of information, one can avoid the controlled- 
phase gates by again using a dual-rail encoding, as in 
Sec. |3.3[ thereby halving the number of qubits that can 
be stored. 

This storage device can be used as the basis for a quan- 
tum computer, where all control is mediated through spin 
1. Such a device is referred to as a Universal Quantum 
Interface [531 H2TH5] . The formulation of the control se- 
quences for this scheme is particularly efficient [JB], al- 
though many other XX model spin chains can also be 

used naim. 



6. Other Fermionic Models 

In this section, we have discussed the properties of mul- 
tiple excitations on a spin chain. This was achieved by 



f A state \ip n ) becomes entangled, via a controlled-phase gate, with 
any state that is still on the chain at the time of removal, and 
which was placed on the chain after \tp n ) was. All controlled 
phase gates that were enacted by previous removals remain un- 
affected. 
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mapping the XX model into a model of non-interacting 
fermions hopping on a nearest-neighbor chain, and prov- 
ing transfer in these systems. However, the class of spin 
models that can be mapped into the same set of non- 
interacting fermions is much broader, which must mean 
that perfect state transfer is achievable in these systems, 
retaining the properties of transfer rate, independence of 
initial state etc. [2"Ul 122] . 

We show this by considering the class of Hamiltonians 



If we define new variables r\kn — {gkn + hkn)/v2 and 
Xkn = {9kn — hk n )v2, then these equations can be suc- 
cinctly expressed as a matrix equation 



li k 



Vk 

Xk 



A + B 
A- B 



Vk 

Xk 



so the combined vector 



Vk 

Xk 



is an eigenvector of a 



N-l 



N 



2N x 2N Hermitian matrix. We can write that 



W = 2 E J n ((l+ln)X n X n+1 +(l- ln )Y n Y n+1 )-J2B n Z n . 



This includes our original XX model as the special case 
7„ = for all n, and the transverse Ising model, consid- 
ered in [201 has j n = 1. As before, the first step is 
to perform the Jordan- Wigner transformation, 



n .JW 



N-l 

E 

71=1 



N 

! E 

n=l 



+2 y ] B n ala n . 



Following [48 , this can be written in a general form of 

9 Bnrn(®n&rn ^m^n) 



H 



general 



N 

E 

n,m—l 



Anm®n®m 



where A is a real, symmetric, matrix (A — A T ) and 
B is real and anti-symmetric (B = —B T ). We want 
to diagonalise this system by introducing new fermionic 
operators ct such that 



N 



(14) 



k=l 



and these new modes obey the canonical commutation 
relations 



} = o 
} = s n 



Our intention is to expand these new modes in terms of 
the old modes, 



E^" " + ^knO-r, 



In order for Eqn. ( 14 1 to hold, it must be that 

^ ^general] = 

which gives a set of simultaneous equations 

m 

l^k^kn ~ ^ ^ Qkm^mn ^km^-rtin- 



E 



VnkVmk 



' Smkhnk 
XnkXmk 



and, similarly, 



k 



QnkQmk ~t~ ^mk^nk 



^ VnkVmk XnkXmk- 
k 

However, these correspond directly to the orthogonality 
relations for the eigenvectors, and hence the canonical 
commutation relations hold. 

Let us take, for instance, the special case of 7 = 1, and 
a nearest-neighbor interaction. This corresponds to the 



transverse Ising model, as considered in 
index the basis elements of 



If we 



M = 



A + B 
A-B 



by \n) for n = 1 . . . 2iV, then we have 

M \ N + n) = B n |n) + 2 J„_i | n - 1) 

M |n) = B n \ N + n) + 2 J n \ N + l + n ) 

This is exactly the nearest-neighbor hopping Hamilto- 
nian of 2N qubits. Hence, we can take any perfect state 
transfer solution for 2N qubits with magnetic field and 
coupling strengths K n and identify B n = K% n _x and 
2J n = K2n, and this is guaranteed to give perfect state 
transfer in this new model. For instance, if we start from 
a state |0), defined such that Ck |0) = for all k, then we 
can create a state a |0) + j3a\ |0). As before the |0) com- 
ponent is an eigenstate, and remains unchanged. The 

1) in M, and gets 

t 

N 



N 



state a\ |0) corresponds to |1) 

transferred to \N_) + \ 2N ) , which therefore outputs ajy |0) , 
so the state is perfectly transferred. 



4. NOISE TOLERANCE 

In a practical situation, no matter how well our system 
is engineered, it will never be perfectly protected from an 
environment (although one might hope that restricting 
access to the chain to just the ends would allow us to 
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protect the bulk of the chain to a large degree). Thus, 
it is worthwhile to make some consideration of errors; all 
this work would be for nothing if the moment a small 
degree of noise is introduced, the quality of transfer were 
to drop radically. 

The presence of static defects, i.e. manufacturing er- 
rors, is a well-studied problem, resulting in the phe- 
nomenon of Anderson localization [49J for sufficiently 
large defects. Small, perturbative, errors lead to per- 
turbative errors in the arrival fidelities [50] . Instead, we 
consider dynamical noise processes, singling out two inde- 
pendent models, local dephasing and local bit flips. The 
individual behaviors of these can be expected to be po- 
tentially very different, although we do not expect either 
to be excessively harmful. In the following subsection, 
we will discuss the case of dephasing noise. The case of 
spin flips is closely related to that of transfer rates and 
the initial state of the system (Sees. 3.4 and 3.3), where 
we learned that the effect of the bit flip is very small 
if the transferred state is not present on that particu- 
lar spin at that particular time. Unfortunately, rigorous 
calculations for bit-flip noise are extremely complicated, 
and not informative, as they necessarily involve the com- 
plete 2 -dimensional Hilbert space of the system. More 
detailed examinations have been made in |511 I52| . 

Ultimately, in the implementation of a quantum com- 
putation, the computational qubits are encoded into er- 
ror correcting codes. If one were to send only one phys- 
ical qubit along the state transfer chain at a time, the 
end result would be independent errors on the physical 
spins, which error correcting codes, and the theory of 
fault-tolerance fo3~|, are well tuned to handling. There 
would be no memory effect since we know how to trans- 
mit through a chain without any effect due to the initial 



evolves to 



state of the chain (Sec. 3.3), which is the only way any 
initially independent noise could potentially build up any 
correlation. 

Of further interest is the fact that it has recently been 
shown that for general spin networks, where perfect state 
transfer is impossible in a coherent way, the transfer abil- 
ity can be vastly enhanced by the presence of noise, and it 
seems that this may have significance for some biological 
processes [54] . 



J t 



where w x denotes the Hamming weight of the bit string 
x, used to enumerate the set of spins on which the Z 
rotations have been applied. The Hamiltonian evolution 
has been incorporated into 

Z x = e- iH ^-^Z x e- iHt . 

If we denote the n th bit of x by x n , then 



z x m =e 



-iHt 



I*) - 2/3 J2 lne- lH{t 



o-t) 



Our aim is to determine how well the state is transferred, 
by calculating the fidelity, 

F = (ip\ Tri..jv_i (pout) \ip) , 

and then averaging over all possible input states \ip) to 
give (F). Ultimately, this yields that 



(F) = l 



2p(2-p) , 2p(l-p) 



N 



Trivial bounds can be applied for the |7 n (0| 4 to give 

2p(2-p) 2p(l-p) 2p 
3 3N - \ ' _ 3 > 

where the upper bound corresponds to the case of simple 
storage without transfer. We note that at any typical 
time during the transfer process, the state will be sig- 
nificantly spread out across the chain, and the fidelity is 
expected to be closer to the bottom end of the range, 
unless one uses a tactic such as encoding the state into 
a wavepacket which remains well localized during the 
transfer [8]. 



2. Independent Baths 



1. Dephasing Noise 

Dephasing noise, the random application of Z errors 
with probability p, independently on each lattice site, 
could be expected to be quite harmful since the relative 
phases that build up between sites during the state trans- 
fer are critical to the perfect transfer. For the sake of an- 
alyticity, we shall restrict to a much simpler model where 
dephasing noise is only applied at one time t (0 < t < to). 
Thus, our initial state 



Another reasonable noise model, discussed in [5T], is 
where each spin of the chain is coupled to a localized set 
of spins by a similar interaction. While we will consider 
the aspects of the proof necessary for the class of chains 
that we are interested in, note that the proof of |51j is re- 
markably general, and essentially any Hamiltonian that 
splits into excitation subspaces will exhibit almost iden- 
tical properties. 

Our aim, as previously, is to transfer states through a 
specific chain, 

N-l N 

Hs = \ J n (X n X n+ i + Y n Y n+ i) — £^ B n Z n , 

n— 1 n—1 
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FIG. 3 The system spins (large), intended for perfect state 
transfer, are coupled to independent baths (small). 



but each of the spins n will also be coupled to an inde- 
pendent bath composed of M n spins as depicted in Fig.[3j 

M„ 

Hb,u = \ ^ 9m( X n X (n,m) + Y n Y {n,m))> 
m— 1 



yielding a total Hamiltonian 



JV 



H — Hs + ^ Hb, ri- 



lt is useful to observe that if an excitation starts on one 
of the system spins n, then while it can hop arbitrarily 
to other system spins, if it appears on a bath spin, it can 
only appear in the states 



G,. 



(n, to) 



1 M n 

22,91 

m—l 

where G\ = J2m=i(9m) 2 ES- This means that a lot of 
the complexity of the bath spins can be eliminated if we 
restrict to the single excitation subspace, with 



where n' is a single effective spin which acts as the bath 
for spin n. We now assume that all the G n are equal 8 . 
The eigenvectors of the system in the absence of the bath 
can be written as 

N 

Wn) = ^m,n \l!±) ■ 

n=l 

allowing us to make an ansatz for the eigenstates of H' = 

H S +J2n H B,n : 

1 N 

|A^> = - 7f Ewi^) + (-i) fe i^))- 



s This makes no assumptions regarding the individual gjjj, so one 
intriguing possibility would be to deliberately introduce one extra 
spin at each site whose g n can be controlled in order to tune the 



These satisfy 

<A^| H> \\ l n ) = S m , n ((-l) k GS kJ + iA n ), 

so although H' is not instantly diagonalized, it is reduced 
to 2 x 2 sub-blocks which are diagonalizable to give the 
eigenvectors and energies. 



E, 



i(A„ + (-l)V4G 2 + A2) 



Having solved for these, we can now turn our attention 
to how well this system implements state transfer. The 
accuracy of transfer of the system in the absence of the 
baths at time t is given by 7jv(*) = 7*(*o — t). In the 
presence of the baths, one can show two results. Firstly, 
in the weak coupling limit, the first correction to 7jy(t) is 
only of order G 2 . In the strong coupling limit, one finds 
that 

Y N (t) = cos(Gi) 7J v(t/2), 

showing that the transfer takes a factor of 2 longer, and 
is modulated by a fast oscillating term. This means that 
extremely high fidelity, or even perfect, state transfer is 
still possible! 



5. PERFECT STATE TRANSFER IN HIGHER SPATIAL 
DIMENSIONS 

With the solutions for perfect state transfer in ID, it is 
straightforward to construct schemes for state transfer in 
larger spatial dimensions. Let us consider two different 
coupling schemes, one for N qubits with nearest-neighbor 
couplings { J n } and magnetic fields B n = 0, and a second 
one for M qubits, nearest-neighbor couplings {K n } and 
magnetic fields B n = 0. Both have the same transfer time 
to- Now consider an N x M lattice of qubits, indexed by 
their positions (x,y). A Hamiltonian of the form 

N-l M 

ff 2D = \ J n( X (i,j) X {i+l,j) + Y {i,j) Y (i+l,j)) 

i=l j=l 
N M-l 

1=1 j=l 

performs state transfer between any diagonally opposite 
points and (N+l — i, M + the two dimensions 

behave independently, and the eigenvectors of the first 
excitation subspace factorize into the form 

N M 
i=l j=l 



where the eigenvectors for the two schemes are 



= V M \ M 



I) 
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The eigenvalues satisfy 



A 



N.M 



Aiy+A 



M 



The Jordan- Wigner transformation does not apply to 
H 2 d, so perfect state transfer is not guaranteed in the 
presence of multiple excitations for the XX model. Thus, 
the question of whether full state mirroring can occur on 
2D lattices still remains open, although if one assumes 
certain symmetry conditions on the coupling strengths, 
it can be shown that mirroring is impossible 55 . On 
the other hand, were one to use a bosonic model, such as 
that of coupled Harmonic oscillators, this transfer must 
work in all excitation subspaces. 

The described technique is not limited to combining 
two ID chains into a 2D Hamiltonian. The introduction 
of the third chain allows the construction to be extended 
to 3D, and so on ad infinitum. The special case of apply- 
ing this protocol with the chain of two qubits generates 
a uniformly coupled hypercube, as studied in [2J, and is 
the basis (or special case) of several of the perfect transfer 
network variants [21 [TS] . 



6. TRANSFORMATIONS AND MANIPULATIONS 




FIG. 4 Two parallel chains, joined at the start, can be used 
for entanglement generation. 

By selecting the J n and B n so that they correspond to a 
perfect state transfer setting, we are guaranteed that the 
Hamiltonian evolves states |n) into \ N + 1 — n ) in time 
to, no matter what these states are. As such, we can 
construct a plethora of different protocols. 

The simplest of these is involved in producing a beam- 
splitter action. We introduce a new system composed of 
two parallel chains (a and b) of spins, joined at one end, 
spin 1, as depicted in Fig. |4j The Hamiltonian of this 
system can be written as 



The motivation that is often given for the study of 
perfect state transfer is that these devices might be real- 
ized in quantum computers, simplifying some tasks. Of 
course, state transfer is not the only task that one might 
want to realize. Having laid the ground work of perfect 
state transfer, however, it is vastly simpler to discover 
how these other 'gadgets' might be realized, as we will 
momentarily examine. 

Nevertheless, the engineering constraints on these 
modulated chains are quite stringent, and may not sim- 
plify anything at all. In that case, is all that we have 
learned wasted? Certainly not! It turns out that the 
study of perfect state transfer is a fantastically useful con- 
structive technique to design all sorts of schemes based 
on Hamiltonian evolution, whether these be constructive 
or destructive effects. Perhaps one wants to show how 
evolution of a Hamiltonian can implement a quantum 
computation, or perhaps one wants to design a pertur- 
bation or interaction with an environment that disrupts 
a coherent effect (thereby giving bounds on the range 
of usefulness of that effect) [S5]. The details of these 
individual calculations will not concern us here, but in 
Sec. 16.41 we will describe the basic method that one uses. 



1. Interferometry and W-state Preparation 

We can view the action of the Hamiltonian Hi (H act- 
ing on the first excitation subspace) as a discrete hopping 
along a chain of N orthogonal states {|n)}; 

Bi |1) + Ji |2) n = 1 

Hi \n) = { J N -i \ N-l) +B N \N) n = N 

J n -i \ n — 1 ) + B n |n) + J n | re + 1 ) otherwise 



Sn=2 Sig{a,6} B n Z n j + J n ( X n iX n +i i + ^ri.i^ri | 1 ,i) 

+B1Z1 + ^ E l e{a M ^ X ^ + 



Using the position basis as before, i.e. |2, a) indicates the 
presence of an excitation on spin 2a and nowhere else, we 
can write a set of orthogonal states to be |l) = |1) and 

I") = -^(Ka) + \n,b)). 

One can then calculate the effect of the new Hamiltonian 
on these states, and see that it maps to the action of Hi. 
As a result, if we start with a state |1), and allow it to 

N), which 



evolve for a time to, we end up with the state 

is a maximally entangled state between two spins. This is 
easily generalized to creating a W-state on M qubits sim- 
ply by introducing M parallel chains, and joining them 
all to spin 1 with the coupling strength Ji/y/M. These 
manipulations were originally expounded in [7J [ST] , 
although this idea has subsequently resurfaced a number 
of times [58j[59]. 

For the special case of N — 2, but arbitrary M, there 
is a neat trick that could be quite helpful. We can readily 
enumerate the other eigenvectors of the system, 



\W k ) 



a 2TTikj/M 



j<£{a...m} 



and further note that the output from the evolution 
(| Wo)) can be converted into any of these via the ap- 
plication of phase gates. This means that once the to 
evolution has completed, application of local unitaries 
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-J 2 /V2/ \j 2 /V2 -\H2f \j 2 /V2 

U - klV2\ _ _ /J,/V2' 







and in the perfect state transfer time, this evolves to 
cos0 |l+) - sin 6» |l-> = cos(26>) |1) + sin(26>) |27V + 1) . 

Thus, by selecting cos(26') = 1/V2, the chain creates the 
maximally entangled state (|1) + \ 2N + l) )/\/2 between 
the furthermost points of the chain [50] . 



FIG. 5 A spin chain network (TV = 6) that acts like an inter- 
ferometer, detecting the phase <f> of the central coupling via 
the probability of which of the output ports an input excita- 
tion arrives at. 

can trap the entangled state on the spins 2. Indeed, this 
is the basis of the routing networks in [T5] . 

In [55] . these ideas were taken further, constructing 
networks whereby unitary rotations can be performed 
while states are transmitted along the spin chains. For 
example, an interferometer can be constructed using the 
network depicted in Fig. [5] 



2. An Alternate Solution for Entanglement Generation 

Let us consider a chain of odd length, 27V + 1, that 
is capable of perfect state transfer, \n) —> | 27V + 2 — n ). 
We can rewrite this chain as two effective chains given 
by states in the symmetric subspace, {|rt + ) = + 
| 27V + 2 - n ))/y/2} and | 7V + 1) , and the antisymmetric 
subspace {|n_) = (|n) - | 27V + 2 - n) )/V2}. The first of 
these chains is described by the interaction 

AT-l 



)(7V + 1 + | + |7V+1 + )(7V + |) 



)<n + l_| + |n + l_)<n_|) 



n=l 
JV+1 

+ ^2 B n \n 

n=l 
+V2J N (\N A 

and the second is 

AT-l 

ff_ = £j n (|n. 

71=1 
N 

+ Y J B n \n-){n_\. 

n=l 

Furthermore, we know what evolutions these must give 
in the perfect transfer time; |n+) — > |n+) and |n_) — > 
— |n_). Now let us consider altering the chain slightly, 
J/v = \/2cos9Jn, Jn + i = V2 sin 9 Jn+i (recall that 
Jn+i = Jn)- We can now apply exactly the same 
analysis as before, except now we have that = 
(cos6\n) + sin6>|27V + 2-n))} and {|n_) = (sia0\n) - 
cos 9 1 27V + 2 — n))}, leaving H± unchanged, and conse- 
quently leaving their action unchanged. Hence, were we 
to start our system in the state |1), this can be repre- 
sented as 

|1) = cos6>|l+) +sin6»|L) 



3. Signal Amplification and GHZ-state Preparation 

Another simple protocol that might be useful in a 
quantum computer is to amplify a signal before measure- 
ment |21j . The idea would be to take an unknown state 
= a |0) + /3 |1) and produce a GHZ state of the form 
a \0)® N + fj 11)®^, yielding a net magnetization that is 
more readily detected by a measuring device. The sensi- 
ble approach here would seem to be to demand that the 
state jO)®^ should be an eigenstate of our Hamiltonian 
f/amp, and that the states 



\ St n 



10) 



0N-n 



should form the steps that we need. One can check that a 
Hamiltonian term K m = ^X m (l — Z m ^iZ m+ i) satisfies 

10)®" = and 



K„ 



m — 1 = n 




(Zn+i is just taken to be 1), so K m+1 replaces the role of 
the term \{X m X m+1 +Y m Y m+1 ). Similarly, Z m Z m+l re- 
places Z m . Thus, we can write down the required Hamil- 
tonian as 



N — 1 

^amp ^ ^ Jn^n 
n=l 



+ 1 



N 

^ ' B n Z„Z, 

n=l 



n+1- 



Fig. ([6| gives a plot of how well the signal gets amplified 
as a function of time using the set of couplings from the 
analytic solution of perfect state transfer. 

It turns out that one can show a unitary transforma- 
tion between H and -ff a mp [21] , not just their first ex- 
citation subspaces. For example, they have an identi- 
cal subspace structure, where the required commutation 
property of H &mp is 



N 

^amp 5 ^ . Z. 
n=l 



n+1 



so the number of boundaries between sets of |0)s and 
|l)s on a chain is preserved. As a result, one can apply 
the results about multiple excitations to this system as 
well. One such result is that a single excitation on a spin 
n > 2 corresponds to a pair of excitations |n — 1, n) on 
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1 2 3 

Time (fraction of state transfer time) 



FIG. 6 Signal amplification for a system of N = 100 qubits 
using J n = \J n(N — n). Perfect amplification occurs at times 
(2n + l)irj and revival of the initial state occurs at times 
2nto- The solid line indicates the probability of the signal 
(initially |l)) being amplified to the desired final state |a^- 

The dashed line shows the average signal strength (as a frac- 
tion of the maximum strength), and the gray line plots the 
probability that more than half of the spins have been nipped. 
Both the dashed and gray lines are largely independent of the 
number of qubits in the system. 

H , and gets perfectly mirrored to spin N + 2 — n. There- 
fore, this chain is capable of performing perfect state 
transfer. Moreover, because we are using the effective 
two-excitation subspace, this behaves like the encoding 
of Sec. 



3.3), which does not require preparation of the 



initial state of the system and yet, in this case, does not 
require any encoding/decoding of the quantum state to 
be transferred. 



4. Universal Quantum Computation 



introduce an ancillary system known as the clock. This 
clock will have N orthogonal states denoted by \n) c . The 
sequence of states that we wish to step through are now 



® (U n - 



so that by initializing the clock in the state |l) c , in time 
to it will arrive at the state \N) C , and the computation 
will have been implemented. It remains to write down a 
suitable Hamiltonian, 



JV-l 



H comp = ^ J n \n + 1) (n\ c ® U n + h.c. 
n=i 

N 

-J2 B n\n) (n| c ®l, 



where h.c. denotes the Hermitian conjugate. Thus, 
Hamiltonian evolution is as powerful as quantum com- 
putation. Such a construction, using perfect transfer, 
was originally introduced by Peres |61| as a modifica- 
tion of Feynman's scheme for Hamiltonian-based quan- 
tum computation |62j . There are a number of ways in 
which this can be extended. For example, one can im- 
pose further structure on the Hamiltonian, such as a ID 
translationally-invariant and S'L r (2)-invariant structure, 
and still recover that arbitrary quantum computations 
can be implemented H5J [531 IM] ■ These constructions can 
also make use of cellular automata for the implementa- 
tion of the computation, which means that the gates to 
be applied are just repetitions of the same basic gate, 
and the computation is entirely determined by the ini- 
tial state of the system. Alternatively, one can use these 
Hamiltonians to understand how difficult it is to find the 
ground state energies d5J I55r[58] . or to discuss the en- 
tanglement properties of ground states [55] , 



Ultimately, we may be interested in how far one can go 
using these simple constructions. Almost by definition, 
a Hamiltonian evolution cannot be more powerful than a 
quantum computer 11 . So, is its power the same, or less? 
For now, we will not be concerned with constructing a 
Hamiltonian which is local with regard to any particu- 
lar spatial distribution; we are merely demonstrating the 
general technique. Evidently, one can write down a quan- 
tum computation as a sequence of unitaries U\ to t/jy— l 
which should be applied to the initial state IV'in)- How- 
ever, it is not necessarily true that the state after n of 
these unitaries has been applied is orthogonal to the state 
at any other time step m, although this is required for the 
Hamiltonian evolution. To circumvent this problem, we 



7. CONCLUSIONS 

In this article we have reviewed to concept of per- 
fect state transfer, and derived the necessary and suffi- 
cient conditions for a nearest-neighbor XX or Heisenberg 
Hamiltonian. In the case of XX Hamiltonians, we have 
explored the properties of the transfer in the presence 
of many excitations which, amazingly, do not hinder the 
effect. We have also demonstrated how a firm grasp of 
perfect state transfer forms the basis of a constructive 
technique with a wide range of applicability, underlin- 
ing the importance of this protocol, even if perfect state 
transfer chains are never realized in a laboratory. 



h The Hamiltonian evolution can be split into a number of small 
steps, and a Trotter decomposition allows each step to be written 
as n-qubit gates for an n-body Hamiltonian, and these can be 
further decomposed into a standard universal gate set within the 
circuit model of Quantum Computation 53 . 
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